TIME-EVOLUTION OF THE POWER SPECTRUM 
OF THE BLACK HOLE X-RAY NOVA XTE J1550-564 
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ABSTRACT 

We have studied the time evolution of the power spectrum 
of XTE J 1550-564, using X-ray luminosity time series data 
obtained by the Rossi X-Ray Timing Explorer satellite. A 
number of important practical fundamental issues arise in 
the analysis of these data, including dealing with time-tagged 
event data, removal of noise from a highly non stationary 
signal, and comparison of different time-frequency distri- 
butions. We present two new methods to understand the 
time frequency variations, and compare them to the dynamic 
power spectrum of Homan et al. [^] All of the approaches 
provide evidence that the QPO frequency varies in a sys- 
tematic way during the time evolution of the signal. 

1. INTRODUCTION 

Recent years have seen dramatic developments in our ability 
to study various parts of the universe in wavelength regions 
inaccessible from the surface of the Earth and using new 
methods of data acquisition, primarily from NASA satel- 
lites. Among the most fascinating are X ray data that are 
suspected to be emanating from black holes. The X-ray lu- 
minosities of such sources are typically highly variable in 
a quasi-random way, and most investigators take the power 
spectruntf\ of the time series data as the diagnostic tool of 

1 Throughout, we loosely use the term spectrum, but this function 
should be carefully distinguished from the energy spectrum, that is the lu- 



choice. In a number of cases the spectra show peaks at 
well defined frequencies, but with relatively broad profiles 
- in contrast to the very narrow peaks of say pulsars, which 
are strictly periodic. These broad power spectral features 
indicate a roughly periodic behavior, called quasi periodic 
oscillations, or QPO's. A number of investigators have re- 
ported time variations in these power spectra. 

Variations in the spectra are of course of immense inter- 
est because they indicate possible changing physical condi- 
tions and hence give basic information about the presumed 
accretion process - the flow of matter from a disk of gas 
surrounding the black hole into its deep gravitational poten- 
tial well. The basic presumption is that the variations are 
connected with irregularities of the flow, which is presum- 
ably turbulent. Progress in this field is crucially dependent 
on proving that the power spectra are indeed evolving, and 
studying the properties of such variation. 

This paper reports on new tools for studying such spec- 
tral evolution, applies them to data on a well known black 
hole candidate, namely XTE J 1550-564, discovered and stud- 
ied using NASA's RXTE Satellite. The next section briefly 
describes the standard tool used by astronomers, with which 
our methods are to be compared. 



minosity as a function of the energy or wavelength of the radiation. 



2. DYNAMIC POWER SPECTRA 

Astronomers usually approach the question of whether the 
harmonic and other frequency content of time series are 
changing over time by computing dynamic power spectra 
(sometimes called spectrograms) from the time series data. 
The point of view is that once a spectral feature has been 
detected and studied by means of ordinary power spectrum 
analysis [g], then time-frequency analysis [[[J is warranted, 
and the dynamic power spectrum is the tool of choice. 

We surveyed the Astronomical Journal and found 15 pa- 
pers using dynamic power spectra. In all but two cases the 
method of computation is essentially not stated, but we pre- 
sume that the authors have computed the ordinary power 
spectrum in a sliding time-window. This cavalier attitude 
- as though the method is so obvious and universal as to 
be barely worthy of mention - implies that astronomers are 
largely unaware of the shortcomings inherent in the dynamic 
power spectrum 

We have developed two methods that we feel are im- 
provements over the usual tools used by astronomers, but 
are nevertheless suited to the data modes in use in astron- 
omy. 

To exercise our methods we have obtained the raw data 
used in generating the "Dynamic Power Spectrum" in Fig- 
ure 19 of [§]. This image shows a QPO initially located 
around 10 Hz, apparently followed by a transition to lower 
frequencies, about 5-6 Hz. The picture is particularly in- 
teresting because it suggests the existence of a time-varying 
frequency with a non trivial behavior. 

We have been able to duplicate this figure. Since the 
original data are photon arrival times, one has first to esti- 
mate the underlying light intensity (rate of photon arrival). 
This was apparently done by binning the data - i.e. by di- 
viding the time axis in fixed intervals of constant length T s 
and counting the photons present in each. We deduce that 
T s was i seconds, and that the instantaneous spectrum was 
estimated with a Welch's periodogram. Figure [y shows our 
reconstructed dynamic power spectrum, which is essentially 
identical to that in [g]. 

3. A TIME-FREQUENCY INVESTIGATION OF 
THE X-RAY DATA 

Although binning is an intuitive and universally accepted 
density estimation procedure, many other methods have more 
appealing properties [H], especially considering that some 
additional signal processing has to be done on the estimated 
density. A famous approach is the kernel method that con- 
sists in averaging the number of events (photons) contained 
in a relatively narrow window that moves continuously along 
the data. The average is weighted with a normalized kernel 
that is more or less concentrated about the center of the win- 



IS 
30 
35 

i 30 




500 L0O0 L 500 3000 3 500 3000 1500 

t(.l 



Figure 1: Reconstruction of the "Dynamic Power Spec- 
trum" presented as Fig. 19 in [Q], The presence of a fre- 
quency component that makes a transition between 10 Hz 
to 5-6 Hz is clearly visible. 

dow. Typical kernels are hence Gaussian, Hanning, triangu- 
lar, etc. The advantages in applying the kernel methods are 
several: 

• no dependence on phasing of the bins relative to the 
signal (the dependence on bin size is replaced with 
that on the width of the kernel) 

• the resulting density inherits properties of the window 
(e.g. it can be differentiated to the same order as the 
window itself) 

• noise (and quantization) are suppressed 

We therefore decided to apply the kernel method, using 
a normalized Hanning window. 

Inspection of either density estimate indicates the pres- 
ence of slowly varying components, on time scales > 0.3 
second (corresponding to frequencies below 3 Hz). To sup- 
press this feature, with the objective to highlight the pres- 
ence of the QPO, we applied a high-pass filter to the data, 
using the forward-backward technique with a simple elliptic 
filter. This method eliminates the nonlinear phase distortion 
introduced by the filter that can produce undesired effects in 
the time-frequency plane. 

At this point, because of the high noise level (due al- 
most entirely to Poisson counting fluctuations), we estimate 
the instantaneous spectrum by applying the so called "slid- 
ing estimator" \^ to the density. This estimator is the nat- 
ural extension to the spectrogram [fit], and is computed by 
evaluating Welch's periodogram of a windowed version of 
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Figure 2: Time-frequency plot of the light intensity emanat- 
ing from the black hole obtained using the sliding estimator 
technique. 

the signal. The window is slid continuously along the sig- 
nal, and hence the dynamic power spectrum can be seen 
as a subcase of this method. Fig. ^ represents the Time- 
Frequency plane obtained with this technique. (The more 
expressive color image can be found on the CD rom of the 
proceedings or at the web page helinet.polito.it/sasgroup/ 
sliding.htm). The joint use of kernel method, highpass filter- 
ing, and sliding estimator confirms the presence of a time- 
varying spectrum in an outstanding way. 

In particular it is possible to notice the first concentra- 
tion of the QPO around 10 Hz and then the transition to a 
central frequency around 6 Hz. A more careful observation 
points out two interesting facts: 

1 . the transition is a "complex" event, because it seems 
to involve a frequency bifurcation around t = 800 
seconds; 

2. the QPO, especially after the transition for t > 1200 s 
is well represented by what is called a one component 
signal - that is, a time-varying frequency with slowly 
varying amplitude. This can also be noticed in the 
first part of the QPO, but with a smaller intensity. 

It is also useful to observe how the noise is reduced, espe- 
cially at higher frequencies, by the application of the kernel 
method, and how the QPO is highlighted by the highpass 
filtering. 



4. A FURTHER DEVELOPMENT IN THE 
ANALYSIS OF THE X-RAY DATA 

In the first part of this paper we have presented a nonpara- 
metric estimation of the instantaneous spectrum of the in- 
tensity of light. This analysis has pointed out that the QPO 
seems to be made by a time-varying frequency that also 
shows a transition from the 10-11 Hz band to the 5-6 Hz 
band. To accomplish this analysis, we adopted on purpose 
a nonparametric approach. That is, we did not make use of 
any model for the data, because still we were not sure of the 
possible presence of oscillations in the signal. 

Now that the result obtained by using the sliding esti- 
mator points out with a good probability the existence of 
this oscillation, we intend to highlight it by doing a further 
signal processing on the estimated time-frequency distribu- 
tion. In doing so we adopt a semi-parametric approach, us- 
ing both additional nonparametric methods and also a para- 
metric model for the background noise that will help us ex- 
tracting the useful information. 

For convenience we call P(t, ui) the estimated instan- 
taneous spectrum obtained in the previous section, that we 
set as the starting point of the new analysis. This analysis 
is based on three main steps: noise whitening, thresholding 
and amplitude normalization of the QPO. A description of 
the steps follows, in which we indicate with Pk(t, u>) the 
time-frequency distribution obtained at the fc-th step apply- 
ing the proposed signal processing. 

4.1. Noise Whitening 

One of the characteristics of the QPO that we have high- 
lighted with the sliding estimator technique, is that it is em- 
bedded in a high level non-white noise, as can be confirmed 
by visual inspection of Fig. ^. This noise does not have 
a flat spectrum, and our aim here is to whiten it, that is to 
force it to have a flat behavior. 

Of course we know that in what we call noise, that is 
the spectral region outside the QPO, there might be other 
important features of the black hole candidate (for example 
other low energy QPOs), but as said before here we want 
to focus our attention on the main QPO. We accomplish the 
whitening in three steps. 

• We first estimate the power spectral density \N(u>)\ 2 
of the noise by integrating with respect to time the 
time-frequency distribution P(t, ui), 

\N(u)\ 2 = J P(t,ui)dt (1) 

This quantity is called frequency marginal in time- 
frequency analysis, and can be seen as an approxima- 
tion of the power spectrum of the signal. 



Then we fit the obtained noise density |7V(a>)| 2 with a 
cubic polynomial to obtain a fitted spectrum | Nf (u>) \ 2 , 
paying attention to force the fitted curve to be grater 
or equal to a small positive parameter 6. We point out 
that the noise density |iV(a;)| 2 is estimated outside 
the QPO region, whereas the fitted density |7V/(w)| 2 
is valid throughout the entire band. 

Finally we whiten the time-frequency distribution P(t, oj) 
by dividing it by the estimated noise 

P(t,w) 
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4.2. Thresholding 



Now that we have whitened the noise, and we expect it to 
show a flatter spectrum, we threshold the obtained distribu- 
tion Pi(t, uj) to extract the QPO. To do this we act in two 
steps 

• We estimate a threshold 8 by inspecting the frequency 
marginal of Pi (t, u>). 

• We then subtract the threshold from the time-frequency 
distribution, setting to zero any negative value the re- 
sulting difference may have 

P 2 (t,Lu) = P 1 (t,uj)-9, P 2 (t,Lu)>0 (3) 

4.3. Amplitude Normalization of the QPO 

As a final step we want to compensate the instantaneous 
amplitude variations of the QPO that can be noticed from 
Fig. |J and that are probably responsible of the changing 
in time of the energy concentration along the QPO instanta- 
neous frequency. To do this we normalize the instantaneous 
spectrum at each time by diving it by its estimated energy at 
that time, that is 

Pt(t,u) 



(4) 



/ P 2 {t,Uj)dtU 

where the integral of P 2 over time is used as an approxi- 
mation to the energy of the instantaneous signal at that time 
(time marginal). 

In Fig. |5]we plot the final instantaneous spectrum P 3 (t,uj) 
where it is possible to notice the effects of the above pro- 
cedure. The analysis done so far better highlights what 
seems to be the time-varying behavior of the QPO central 
frequency. Our analysis to date seems to elucidate QPO 
variations better than conventional dynamic power spectra, 
but further studies are needed, mainly of the algorithm's re- 
liability and immunity to artifacts. 




1500 2000 
Time [s] 



Figure 3: Extraction of the QPO with a more advanced and 
quasi-parametric signal processing. 
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